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Abstract: We study the dipole picture of high-energy virtual-photon-proton scattering. 
It is shown that different choices for the energy variable in the dipole cross section used 
in the literature are not related to each other by simple arguments equating the typical 
dipole size and the inverse photon virtuality, contrary to what is often stated. We argue 
that the good quality of fits to structure functions that use Bjorken-x as the energy variable 
- which is strictly speaking not justified in the dipole picture - can instead be understood 
as a consequence of the sign of scaling violations that occur for increasing Q 2 at fixed small 
x. We show that the dipole formula for massless quarks has the structure of a convolution. 
From this we obtain derivative relations between the structure function at large and 
small Q 2 and the dipole-proton cross section at small and large dipole size r, respectively. 
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1. Introduction 

The colour dipole model [1,2,3] provides a successful description of deep inelastic scattering 
(DIS) processes in a wide range of the kinematic variables. We consider here electron- and 
positron-proton scattering, 

e ± +p^e ± + X, (1.1) 

at not too high momentum transfers squared Q 2 , Q 2 < 1000 GeV 2 say, and at high energies. 
There, only the exchange of a virtual photon 7* between the leptons and the hadrons has 
to be taken into account. Thus, we study in essence the absorption of a high-energy virtual 
photon on a proton, 

j*+p^X. (1.2) 

The structure functions for this inclusive DIS process were extensively measured at HERA 
[4, 5, 6, 7, 8, 9]. A considerable part of the literature concerning the proton structure 
functions in DIS uses the colour dipole picture as an essential input. In particular, the 
interpretation of the structure function F2 at small Bjorken-x relies heavily on the dipole 
picture. The corresponding results are frequently used in the context of various scatter- 
ing processes in proton-proton collisions at the LHC and in the description of the initial 
conditions of the creation of the quark-gluon plasma in heavy-ion collisions. Also calcu- 
lations for processes at a future lepton-hadron collider [10] are often based on the dipole 
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picture. Obviously, a good understanding of the dipole picture and its consequences is 
very important for all these applications. In the present study we shall consider an aspect 
of the dipole picture that has - in our opinion - not yet received proper attention in the 
literature, namely the correct choice of energy variable in the dipole-proton cross section. 

The idea that the high energy photon in the reaction (1.2) acts in some way like a 
hadron goes back a long time. It has been used for instance since the 1960s in vector 
dominance models, see [11, 12] for reviews. Today, the dipole picture is frequently used 
in analyses of DIS structure functions. There, the reaction (1.2) is viewed as a two-step 
process. In the first step the photon splits into a quark-antiquark pair which represents 
the colour dipole. Subsequently, that pair scatters on the proton, this second step being 
a purely hadronic reaction. For reviews we refer the reader to [13, 14]. In [15, 16] the 
foundations of this dipole picture were examined in detail. The precise assumptions which 
have to be made in order to arrive at it were spelled out. In [17, 18] it was shown that 
already the general formulae of the standard dipole approach allow one to derive stringent 
bounds on various ratios of structure functions. These bounds were used to determine the 
kinematic region where the dipole picture is possibly applicable. In particular, it was found 
that for 7*p c. m. energies W in the range 60 to 240 GeV the standard dipole picture fails 
to be compatible with the HERA data for Q 2 larger than about 100 to 200 GeV 2 , see Fig. 9 
of [18]. 

The derivations of some of these bounds rely on the dipole-proton cross sections fr( q \ 
where q denotes the quark flavour, being independent of Q 2 . In [16] it has been stressed 
that this Q 2 -independence of the dipole-proton cross section is in fact natural. The correct 
energy variable is given exclusively by W and the functional dependence of should be 

= &to)(r,W), (1.3) 

where r is the transverse size of the dipole. This excludes in particular the choice of 
Bjorken-x instead of W, since this would introduce a dependence on Q 2 in addition to W, 
see (2.2) below. In the derivation [15, 16] of the dipole picture the dipole cross section 
arises from a T-matrix element for the scattering of a dipole state on the proton. The key 
feature of these dipole states is that they consist of a quark and an antiquark described by 
asymptotic states. The dipole states are then independent of Q 2 in the high energy limit. 
Upon a smearing in the relative transverse vector r between quark and antiquark and in the 
longitudinal momentum fraction a of the photon carried by the quark the dipole states can 
be viewed as hadron analogues, whose normalisation is independent of continuous internal 
degrees of freedom. But also the mean squared invariant mass of such smeared dipole 
states is independent of Q 2 at large W . Since any physical cross section can depend only 
on variables defined by the incoming states, the dipole-proton cross section hence cannot 
depend on Q 2 . Nevertheless, the energy variable x - and therefore a Q 2 -dependence - 
is frequently used in popular models for the dipole cross section, such as [19]. Further 
examples for x-dependent dipole cross sections are [20], [21], and [22]. Sometimes also 
other dependencies on Q 2 are introduced by modifying the photon wave functions [23, 24]. 
Other models introduce an impact parameter dependence of the dipole cross section, see 
[25, 26]. For an overview of the physical motivations for various models see for instance 
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[27]. Furthermore, the perturbative gluon density g(x,Q 2 ) naturally depends on x and Q 2 
and its often assumed proportionality to the dipole cross section suggests that the latter is 
also Q 2 -dependent. A subtle point in such a comparison is given by the fact that different 
limits are used for the dipole picture and the double leading logarithmic approximation, 
respectively, see the discussion in [16]. Only very few models for the dipole cross section 
have been constructed that use the correct functional dependence (1.3), among them are 
[24], [28], and [29]. 

Through the interplay of photon wave function and dipole cross section a typical trans- 
verse dipole size is generated, which depends on Q. In this paper we shall investigate 
whether such an effective dipole size can be used to relate different choices of energy vari- 
ables in the dipole cross section. 

We consider then the dipole formulae for the case of massless quarks and show that 
these formulae can be understood as a convolution. This is used to analyse the relation 
of the structure functions and the dipole cross section. For high and for low Q 2 we find 
simple but somewhat surprising relations. 

Our paper is organised as follows. Section 2 reviews the relevant formulae of the 
dipole picture. In section 3 we discuss typical dipole sizes and investigate whether different 
choices for the energy dependence of the dipole cross section may be related to each other by 
effective scale arguments. In section 4 we rewrite the dipole formula for massless quarks as 
a convolution. We present considerations suggesting that the success of using Bjorken-x as 
the energy variable in the dipole cross section is due to the specific form of scaling violations 
in the structure function F<i- We derive asymptotic relations for the general dipole picture 
in the regimes of large and small Q 2 , respectively. Our conclusions are drawn in section 
5. Two appendices contain results used in the main text: In appendix A we derive the 
asymptotic behaviour of the integrated photon densities. Appendix B explains the steps for 
obtaining a simplified version of a VF-dependent model for the dipole-proton cross section 
from the original model [29]. 



2. The dipole picture 

We use the standard formulae for the kinematics and for the definitions of structure 
functions of the reaction (1.1), see for instance [30]. As discussed above, we consider 
Q 2 < 1000 GeV 2 so that it is sufficient to take into account the exchange of a photon. 
Thus, we shall study in the following the absorption of a virtual photon 7* on the proton, 

l*(q)+p(p)^X(p'). (2.1) 

Here the 4-momenta are indicated in brackets. The c. m. energy for this reaction is denoted 
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Figure 1: Basic diagram for the description of the cross sections <jt,l of j*p scattering in the 
standard dipole approach. 



by W, the virtuality of 7* by Q . For these and the other usual variables we have 

r2 _ 1 „\2 

Q 1 



W 2 = (p + q) 2 , 

l2 _ „2 



p-q 
v = , 

2m p 

X = ^L = 91 (22) 

2m p u W 2 + Q 2 -m 2 ' y ' ' 

The proton in (2.1) is supposed to be unpolarised, while the virtual photon can have trans- 
verse or longitudinal polarisation. The corresponding total cross sections are ar(W, Q 2 ) 
and crL(W,Q 2 ), respectively. The F 2 structure function is, with Hand's convention [31] for 
the 7* flux factor, 

O 2 ( m 2 \ 

F 2 (W,Q 2 ) = ^—[a T (W,Q 2 ) + a L (W,Q 2 )](l-x) + oi^L\ . (2.3) 

For small Bjorken-x, x <C 1, this simplifies to 

F 2 (W, Q 2 ) = [a T (W, Q 2 ) + a L (W, Q 2 )] . (2.4) 

In the following we shall use this simpler relation since we shall only consider the region 
x < 1. 

In order to obtain the standard dipole model for the cross sections cjt,l we can relate 
them first to the imaginary part of the j*p —> 7*p forward scattering amplitude. The latter 
is represented as the initial 7* splitting into a qq pair, this pair scattering on the proton 
and the qq subsequently fusing into the final state 7* , see Figure 1 . With the assumptions 
spelled out in detail in section 6 of [16] the diagram of Figure 1 gives in the high energy 
limit 

vt,l(W,Q 2 ) = Y1 J d 2 rw^ L (r,Q 2 )a^(r,W). (2.5) 
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Here the integrated 'photon densities' are 



</S + (a,r,Q) 



A, A' J ° 
A, A' J ° 



(2.6) 
(2.7) 



We recall that a is the longitudinal momentum fraction of the 7* carried by the quark, r 
is the vector in transverse position space from the antiquark to the quark, r = \r\, and A 
and A' are the helicities of q and q, respectively. The total cross section for the scattering 
of the qq pair on the proton is denoted by fr( q \ the 7* wave functions for transversely and 
longitudinally polarised 7* by ip^x^ an< ^ ^AA'^' respectively. A sum over all contributing 
quark flavours q is to be performed in (2.5). Inserting in (2.6) and (2.7) the photon wave 
functions in leading order of a s and assuming the longitudinal momenta of quark and 
antiquark to be much larger than their mass and transverse momenta we get the standard 
expressions 



£ |V>S + («, r,Q) * = a em Q 2 q { [a 2 + (1 - a) 2 ] e^e.r)] 2 + m 2 q [K (e q r)] 2 } , (2.8) 



N, 



A, A' 



£Ug,V,r,Q) 



A, A' 



2N r 



71" 



a cm Q 2 q Q 2 [a(l - a)] 2 [K (e q r)] z . 



(2.9) 



Here N c = 3 is the number of colours, Q q is the charge of the quark in units of the proton 
charge, -Kq,i are modified Bessel functions, and 



e q = \J a{l - a)Q 2 + m\ . 



(2.10) 



For massless quarks q the integrated photon densities w^ L of (2.6) and (2.7) can be cal- 
culated analytically. The corresponding formulae are given in appendix A. 

In (2.5) we have written ot and ol as functions of W and Q 2 , and the dipole cross 
sections as functions of r and W. It is one purpose of this paper to examine the 
consequences of choosing energy variables other than W in the dipole formulae. Choosing, 
for instance, Bjorken-x as variable we get again the formulae (2.5) but with W everywhere 
replaced by x. On the l.h.s. of (2.5) this means, of course, just another and equivalent 
pair of variables. But for the r.h.s. this has drastic consequences as we shall show in the 
following. 



3. Energy dependence of the dipole cross section 
3.1 Typical dipole sizes 

In [16] it was shown that in a formulation starting from the functional integral the dipole 
cross section comes out as independent of Q 2 and depends on r and the energy W 
only. We already noted that in contrast to this, prominent dipole models discussed in the 
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literature introduce a dependence of on Q 2 , typically through the Bjorken-x variable. 
Often it is argued that in the dipole model one has a relation of the kind 

r = |, (3-D 

with a constant C, corresponding to a 'typical dipole size' or 'typical scale'. That is, the 
most relevant dipole sizes r are determined by the scale Q. The reasoning behind this is 
the fact that due to the interplay of the Q 2 -dependence of the photon wave function and 
the r-dependence of the dipole cross section a typical size is generated. Neglecting the 
W-dependence, the masses and all non-perturbative intrinsic scales, the typical size must 
be given by (3.1) for dimensional reasons. Taking (3.1) seriously one might be tempted to 
replace freely in a(r, W) or related quantities r by Q 2 dependencies and vice versa. In this 
section we shall show that such replacements are far from harmless and, in fact, are not 
admissible. 

Let us first see how the relation (3.1) for typical dipole sizes arises. Inserting (2.5) into 
(2.4) we can write the dipole model expression for the structure function F 2 as 

f°° 1 

F2(W,Q 2 ) = QJ dr^2h(Qr,m q r)-^Q 2 q a (<1 \r,W). (3.2) 

Here we define 

h(Qr,m q r) = Qr \ 2 \wfi(r, Q 2 )+w^\r, Q 2 )j . (3.3) 

Z7TQ! em t^q L J 

The dependence of the dimensionless functions Qr 3 w^ L / (a cm Q 2 ) on r, as resulting from 
(2.6)-(2.9), is shown in Figure 2 for the case of massless quarks. The function h is shown in 
Figure 3. For massless quarks, m q = 0, h is only a function of the dimensionless variable 

z = Qr. (3.4) 

The function h(z, 0) has its maximum at 

z = 2.4010. (3.5) 
The behaviour of h(z, 0) for small and large z is 

h(z, 0) oc z for z — > , 

1 (3-6) 
h(z, 0) oc - for z — > 00 . 

In many proposed models colour transparency at small r is implemented by assuming 
<j(9)(r) oc r 2 for r — > 0. For larger r, on the other hand, the dipole cross sections are 
certainly not expected to grow faster than r 2 . Thus Q 2 a^ q \r)/r 2 should be a rather 
smooth function of r. Let us now, for the sake of argument, consider massless quarks only. 
Then we get from (3.2) 

F 2 (W,Q 2 ) = Q J™ drh(Qr,0)f (r,-^j , (3.7) 
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Figure 2: Dependence of the dimensionless quantity Qr 3 w^ L /(a cm Q 2 ) on r for massless (q = 
u, d, s) quarks for three different values Q 2 = 1, 10, 100 GeV 2 ; on the left for transversely polarised 
photons, on the right for longitudinally polarised photons. 
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Figure 3: The function h(z,0) dehncd in (3.3). 



where we define the function / by 



(3.8) 



In (3.7) a smooth function of r, f(r, l/(r 2 W 2 )), is integrated with h(Qr,0) having a maxi- 
mum at r max = zq/Q. It is now tempting to replace r in the smooth function f(r, l/(r 2 W 2 )) 
by r max . In this way we get a modified F2, 



Fr d '(W, Q 2 ) = Q J™ dr h(Qr, 0) / (§, ^) 



(3.9) 
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where we used Q 2 /W 2 « x for W 2 ^> Q 2 . By this trick the effective dipole cross section 
got a Q 2 -dependence. 

Of course, this is too simplistic and does not work since the integral in (3.9) diverges 
at large r due to (3.6). But we can modify the above argument slightly and replace r by 
r max = zq/Q only where r is associated with the energy scale W. In this way we get 

Fr d (W,Q 2 ) = Qj™drh(Qr,0)f(r,^j . (3.10) 

Now the integral is in general convergent and we obtain a dipole formula with a dipole- 
proton cross section depending on r and x. That is, we made the replacement of scales 
(3.1) 

52Q 2 q *M(r,W) — ► r 2 f(r^) = a™ d (r,x) (3.11) 
i 

with C = zq = 2.40. But also other values of C can be envisaged. 

We now show that also the replacement (3.11) is far from harmless and, in fact, modifies 
the structure function in an essential way. Basically this is due to the fact that the function 
h(z, 0) is very broad and is not well approximated by a delta function at zq. 

3.2 Substitution of scales via typical dipole sizes: examples 

We first consider an example where a Q 2 -dependence is introduced via the 'typical dipole 
size' into a dipole cross section that originally depends only on r and W but not on Q 2 . 
Actually there are only few examples for dipole models of the latter kind in the literature. 
We choose here a slightly simplified version of a model proposed by Donnachie and Dosch 
in [29]. This DD model is based on Regge theory and includes exchanges of a soft and a 
hard pomeron. The intercepts of the soft and hard pomeron trajectories are denoted by 
(1 + e s ) and (1 + e/J, respectively. The dipole cross section that we want to consider is 



(3.12) 



The parameter values are 



e/l = 0.42, A = 57.4mb/fm, i? c = 0.22fm, 

e s = 0.08, W = 20GeV, a = 0.346fm. (3.13) 

We take into account only light quarks (it, d, s) and set their masses to zero for simplicity. 
In appendix B we show how this simplified model is obtained from the original, more 
general approach of [29]. Depending on the kinematic parameters the simplifications lead 
to non-negligible deviations from the original model. But the simplified version is sufficient 
for the sake of our argument even if it describes the data only moderately well in some 
kinematic region. We should point out that both the simplified and the original model 
exhibit a rather strong rise of F2 at very small x and large Q 2 that originates from the 
assumed high intercept of the hard pomeron in this model. However, substantial deviations 
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from the Golec-Biernat-Wiisthoff model (see below) occur only in regions in which no data 
are available. 

For the DD model we hence have 

a DD (r, W) = Y, Q*Al(r, W) = r 2 f DD (r, -^j , (3.14) 



where the sum is over the light flavours only and 



/dd r, 



1 



r 2 W 2 



q 

x 9{R c -i 



exp 



( 3.1a) 



( r 2 W 2 Y 



+ 6{r - R c 



r 2 W 2 



(3.15) 



Making the replacement of scales (3.11) in (3.15) gives 



a^ d (r,x) = r'f DD [r, 



a 2 ") 



(3.16) 



with C = zq = 2.40. As discussed above, this replacement changes the cross section 
only via terms which are sensitive to the external scale W, that is the terms in brackets 
with exponents and e s in this case. Figure 4 shows the effect of this substitution on 
the structure function F2 obtained from these dipole cross sections (calculated only with 
light flavours). The solid curves show F2 obtained using the (original) VF-dependent dipole 
cross section <7dd(^ W), while the dashed curves show F2 calculated using the x-dependent 
o"]3j) d (r, x). The results clearly deviate from each other. In particular, the x- and (In- 
dependences are heavily altered. This is also the case when one uses in the replacement of 
scales (3.1) not C = zq but C = 2zq = 4.80, corresponding to larger 'typical dipole sizes'. 
F2 obtained with this latter replacement is shown as the dotted curves in Figure 4. 

Let us now consider the reverse substitution of scales. That is, let us start from an 
x-dependent dipole cross section a(r, x) and make the reverse replacement of (3.11), namely 



a r, x 



- r 2 f (r, r 2 f (r, ee a^(r, W) (3.17) 



with C = zq = 2.40. As an example we consider the dipole model proposed by Golec- 
Biernat and Wiisthoff [19]. This GBW model describes the F 2 data from HERA quite 
well. The dipole cross section of this model is given by 



a GB w(r,x) = ^Q^GBw( r ' x ) = ^2Q 2 q a o 
1 q 

Accordingly, we have 

/gbw (r, -§2) =^Y1 Q l a ° 



1 — exp 



2R (x) 



(3.18) 



1 — exp 



2R (x) 



(3.19) 
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Figure 4: obtained from the simplified DD model for a and modified versions thereof, with 
energy dependencies via 'typical dipole size' substitutions, for three values of Q 2 — 2, 10, 100 GeV 2 . 
Bjorken-x is varied and Q 2 is kept fixed for each plot. Shown are the DD model with original 
W-dependence (solid curves) and modifications of type (3.11) with C = z — 2.40 ('modi', dashed 
curves) and C = 2z = 4.80 ( l mod2', dotted curves). 



For simplicity we consider also here only light (u, d, s) quarks, neglect their masses, and 
choose the following parameter set of [19]: 

/ \ A/2 

R (x)= ( — ) GeV" 1 , fj = 23mb, A = 0.29 , x = 3 • 10" 4 . (3.20) 

The replacement (3.17) gives us a modified dipole cross section OqbwC 7- ' ^0- Figure 5 shows 
the effect of this replacement on Fi- The curves for F% of the original GBW model (solid 
lines) are significantly modified by the substitution (3.17) (dashed lines), in particular 
the dependence on Q 2 and x is altered. Again, this is also the case when using in the 
replacement of scales (3.1) not C = z$ but C = 2zq = 4.80, corresponding to larger 'typical 
dipole sizes', see the dotted lines in Figure 5. 

In this section we have shown on two examples that the substitution of scales r <H> C/Q 
(3.1) in the dipole cross section a, done in either way, alters the structure function F2 
significantly. As could have been expected, the direction of the alteration is opposite in 
the two directions of performing the substitution. (Clearly, if one applies the substitution 
and its reverse subsequently on the same dipole cross section a, the individual alterations 
have to cancel. The effect persists if one applies the two ways of substituting to different 
models which individually describe the data well.) The size of the alteration due to the 
substitution of scales is significant both quantitatively and qualitatively, see Figures 4 and 
5. 

Actually, we can infer already from the results of [17, 18] that the choice of the energy 
variable in a is crucial at least at high Q 2 . There it was shown that any dipole cross 
section of the form a(r, W) fails to describe the HERA data for Q 2 > 100 to 200 GeV 2 . In 
contrast, the GBW model with <5"GBw( r ) a; ) provides a good fit to the data also at higher 
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Figure 5: F2 obtained from the GBW model for a and modified versions thereof, with energy de- 
pendencies via 'typical dipole size' substitutions, for three values of Q 2 — 2, 10, 100 GeV 2 . Bjorken-x 
is varied and Q 2 is kept fixed for each plot. Shown are the GBW model with original x-dependence 
(solid curves) and modifications of type (3.17) with C = zq = 2.40 ('modi', dashed curves) and 
C = 2z = 4.80 ('mod2', dotted curves). 



Q 2 . From this we see that a dependence of a on Q 2 in addition to W can certainly not be 
eliminated or introduced by an effective scale argument of the type (3.1) in the regime of 
high Q 2 without drastic consequences. Let us, indeed, compare ratios of F2 for different 
values of Q 2 for the models, original and modified, discussed above. An illustration of such 
ratios 1 is given in Figure 6 for the example W = 150 GeV, where in addition the general 
bound (10) of [17] is shown, which is valid for any a(r, W). We see that, as expected, 
^GB\v( r '^) as wen as &~DD(r,W) respect the general bounds. In contrast, the bounds are 
violated for both x-dependent cross sections ffGBw( r !^) and <5"gg d (r, x). 



4. The dipole formula as a convolution 

In this section we consider again, for simplicity, only the light quarks u, d, s and neglect 
their masses. Our starting point is the relation (3.7) with / inserted from (3.8). But now 
we leave the choice of energy variable open, 

POO 1 

F2(C,Q 2 ) = QJ o drh(Qr,0)^a(r,0, &(r,£) = 5>2* (ff) ( r >0 ■ (4.1) 

We recall that the function h(Qr, 0) is obtained directly from the photon wave functions, 
see (3.3). Here we set £ = W if we want to consider a VF-dependent dipole cross section 
and £ = x for a x-dependent one. Correspondingly, we consider F2 as function of £ and 
Q 2 - 

1 Note that the curves are calculated from the models and shown here for a kinematic range that is 
slightly larger than that in which data from HERA are available. Corresponding curves restricted to the 
actual kinematic range where HERA data exist can be found in [17, 18]. 
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Figure 6: Ratio F2(W, Qf)/.F2(W, Q%) for two dipole models and modified versions with altered 
energy dependencies via 'typical dipole size' substitutions, for W = 1 50 GeV. Here, Q\ is varied 
and Q\ = 10 GeV 2 is kept fixed. Left plot for the (simplified) VF-dependent Donnachie-Dosch 
model (solid curve) and its modification (3.11) with C — zq — 2.40 (dashed curve). Right plot 
for the original Golcc-Bicrnat-Wiisthoff model (solid curve) and its modification using (3.17) with 
C — zq — 2.4 (dashed curve). Quark masses are set to zero in all cases. The shaded area is excluded 
by the bound (10) of [17] for any Q 2 -m.dependent dipole cross section a^ q \r, W). 



Now we show that (4.1) is a convolution formula. For this we set, as in (3.4), Qr = z 
and define (recall that zq is defined as the position of the maximum of h(z,0)) 

t = ]n(z/z ) , (4.2) 

k(t) = z h(z e T ,0). (4.3) 

We have < z < oo and — oo < r < oo. Corresponding to the behaviour of h(z,0) for 
small and large z, see (3.6), we have 

k(t) oc exp(±r) for r — > =Foo . (4-4) 

The function k(t) is shown in Figure 7. The half width of this function is At = 2.22. 
Now we choose Qq and ro such that 

Qom = zq (4.5) 
and define new variables, replacing Q and r, 

t = HQ/Qo) (4.6) 

and 

t' = -ln(r/r ). (4.7) 
Furthermore, let us define a rescaled and reparametrised 'dipole cross section' 



rro 



r=ro exp(— t') 



(4.8) 
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Figure 7: The function k(t) defined in (4.3). 
The dipole formula for massless quarks (4.1) has now the form of a convolution 

/oo 
dt'K(t-t')S(C,t'). (4.9) 
-oo 

For dipole cross sections which are reasonably behaved for small and large r the integral 
in (4.9) is convergent. Indeed, we find 

S(f , oc exp(-i') for t 1 -> oo, if &(r, f ) oc r 2 for r -> (4.10) 

and 

S(f,0 oc exp[(l - a)t'] for t' ^ -oo, if <r(r, £) oc r a for r -> oo . (4.11) 

In the following we shall assume o < 1. That is, we assume that <r(r, £) increases more 
slowly than linearly with r at large r. In many models the dipole cross section saturates for 
large r which corresponds to a = in (4.11). With these assumptions the function S(£,t') 
decreases exponentially for t' — > ±oo and so does n{t — t'). Thus, the convolution formula 
(4.9) is well convergent. 

We assume now, for the sake of the argument, that for some given £ the dipole formula 
(4.1) is valid for all Q 2 , that is, (4.9) holds for all t. This allows us to draw some general 
conclusions: 

(i) The convolution formula (4.9) by itself places really no restrictions on i*2- We can, 
in principle, invert the formula, using well known techniques from the theory of 
integral transforms, and obtain S(£, t') and the dipole cross section from i 7 ^- However, 
inversion is in general difficult and may constitute an ill-posed problem depending 
on the detailed structure of the functions involved. The physical restrictions enter in 
(4.9) from the requirement of a non-negative dipole cross section, that is, from 

S(£,t')>0 for alii'. (4.12) 
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Suppose now that we have indeed inverted the convolution (4.9) and found a non- 
negative S(£,t'). Then the dipole-proton cross section a(r, £) will appear as an in- 
tegral over Q 2 of F 2 (^,Q 2 ), respectively <tt(£, Q 2 ) + ctl(£,,Q 2 ), using (2.4). Setting 
now £ = x we would find that &(r, x) is related to the j*p cross sections at the same 
x and all Q 2 . This means that the dipole cross section at some fixed r and x would 
appear as (inverse) convolution of the structure function F 2 (x,Q 2 ) at the same fixed 
x but varying Q 2 corresponding to very different hadronic final states X in (2.1). 
Taken literally, the final states would range from the single proton (Q 2 = 0) to multi- 
particle final states with arbitrarily high invariant mass for Q 2 — > oo. This looks not 
so plausible to us from the physics point of view. On the other hand, setting £ = W, 
only 7*p cross sections for final states with the same c. m. energy will appear in the 
integral for a(r, W). This appears much more plausible to us. 

(ii) The function n(t — t') in (4.9) is rather broad and decreases as exp(— t) for t —> oo, 
see (4.4) and Figure 7. Superposing n(t — t') with non-negative weights S(£, if) will 
lead to an even broader function of t on the l.h.s. of (4.9). That is, the decrease of 
-^2(£) Qo e2t ) ex P( — t) can then not be faster than exp(— t) for t — > oo. This means 
that the dipole formula will be in difficulties if, for large Q 2 , F 2 (^, Q 2 ) decreases with 
Q 2 at fixed £ but will work if i*2(£)Q 2 ) increases with Q 2 at fixed £. From this 
simple observation we can already see why, phenomenologically, dipole models with 
the choice £ = x work up to higher Q 2 than models where the choice is £ = W . For 
fixed small x the structure function F 2 {x, Q ) increases with Q 2 in the HERA regime, 
see left panel in Figure 8. For fixed W, in contrast, F 2 {W, Q 2 ) first increases, but 
eventually decreases with increasing Q 2 , see right panel in Figure 8. 
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Figure 8: Q 2 -dependence of the structure function F2. On the left for fixed x = 5 • 1CP 3 , on the 
right for fixed W = 150 GeV. F2 shown here is obtained from the parametrisation of [32], the band 
indicates a la- uncertainty. 

We shall now illustrate the convolution formula (4.9) for the example of the GBW 
model where we have, of course, the choice £ = x, see (3.18) and (3.20). Furthermore we 
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set 

Qo = GeV , r = z /Q = 0.24 GeV" 1 = 0.047 fm . (4.13) 
In Figure 9 we show for x = 10 -3 the original function Sgbw> suitably normalised, 

0.15 



/f'QO 
«)gbw(z,*) := SgbwOM) / drfi(r) = — ffGBW^a;) 



r=ro exp(— t) 



(4.14) 



as well as the result of the convolution, 

/oo 
dt' n(t - t')S G Bw(x,t') = F 2G bw(x,Q 
-oo 



2n Qo 



Q 



=Q exp(t) 



. (4.15) 



The effect of the broadening due to the convolution is clearly visible in Figure 9. We also 
see that within a factor of 2 we have very roughly 



(5 * «)gbw(M) ~ (S J «)gbw(M) for - 3<t<l. 



(4.16) 



This relation actually holds for a wide range in x, with the t-range varying somewhat with 
x. Translating this back to r and Q we find that within a factor of 2 we have 



Q 



0.15 

rro 



r=zo/Q 



for 0.2 GeV < Q < 8 GeV . (4.17) 




Figure 9: The functions (4.14) (red, dashed) and (4.15) (blue, solid) in the GBW model for 
x = 1(T 3 . 



But for large and for very small Q 2 a relation of the form (4.16) respectively (4.17) 
does not hold at all. In the left panel of Figure 10 we show the functions (4.14) and 
(4.15) of Figure 9 but multiplied by exp(i). This gives a clear indication of what happens 
for large t, that is, large Q 2 and small r. We see from the Figure that, for t > — 1, 



- 15 - 



F*2 gbw = e\S * k)gbw and e\S J k) gbw = 0.15 r 2 a GB w\r=r exp(-t) have a completely 
different behaviour. F2GBW rises linearly with t, while ogbw/^ 2 goes to a constant for 
t — > 00, that is r — > 0. The latter corresponds to the assumption of colour transparency 
in the GBW model. Note that the derivative of a linear function of t gives a constant. 
We shall show below that quite generally we expect indeed a derivative relation between 
F 2 {£,Q 2 ) and r" 2 a(r,0 for large Q 2 . 
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Figure 10: The two functions of Figure 9, but now multiplied by exp(t) (left panel) and exp(— t) 
(right panel). 



In the right panel of Figure 10 we show again the functions (4.14) and (4.15) of Figure 
9 but now multiplied by exp(— t). In this way we see clearly what happens for t — > —00, 
that is, for very small Q 2 and large r. For t — > —00 the function e~ l S* k = e _2 *i*2 behaves 
as (—t) whereas e~ l S f k is constant, suggesting again a derivative relation between them. 
Indeed, we shall show below that on general grounds we expect for t —> —00 such a 
derivative relation between gt + cl and a to hold. 



4.1 The large-Q 2 regime 

In this subsection we give a general discussion of the convolution formula (4.9) for large t 
corresponding to large Q 2 and small r. We rewrite (4.9) in the form 



/oo 
dt' K (i - t') exp(t - t')e* ' S(£, 
-00 



where 



e t 'S(d,t')= ^Hr,0 



r=ro exp(— t') 



(4.18) 



(4.19) 



The function k(t) exp(r) is shown in the left panel of Figure 11. Very roughly it resem- 
bles a step function. Suppose now that cr{r^)/r 2 is slowly varying for r — > corresponding 
to exp(t')S(£,t') being slowly varying for t' — > 00. In (4.18) this function is folded with 
n(t — t')exp(t — t') which resembles a step function. The folding produces, therefore, in 
essence the primitive of exp(t)£(£, t). That is, we expect for large t the following relation 
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Figure 11: The functions «;(r)exp(T) and K(r)exp(— r) 



Cooe'S^t)^-^ F 2 (£, Qle 2t ) (t » 1) , 



(Q 2 > 74GeV), 



with 



coo := 2vt(k(t) exp(r))| r ^ 00 = ^ = 0.065. 



7T" 



(4.20) 



(4.21) 



(4.22) 



This is illustrated for the GBW model in the left panel of Figure 12 for the choice x = 10 3 . 
Here the relation (4.20) respectively (4.21) is indeed well confirmed for t > corresponding 
to Q 2 > 10 GeV 2 . 
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Figure 12: Numerical test of the derivative relations (4.20) respectively (4.21) (left panel) and 
(4.26) respectively (4.27) (right panel) for the GBW model for x = 10~ 3 . 



From (4.21) we see that the dipole formula will work well phenomenologically at high 
Q 2 if we choose an energy variable £ such that i^Cd Q 2 ) increases with Q 2 at fixed £ since 
then the cross section a can be positive. This is the case for the choice £ = x for small x 
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but, of course, many other choices of £ would also have this property. On the other hand, 
F2CW, Q 2 ) decreases for fixed W at large enough Q 2 , see figure 8. Therefore, it is clear that 
dipole models with the choice a(r, W) will only be able to describe the structure function 
F^iWi Q 2 ) over a more limited range of Q 2 values. But this limitation at high Q 2 is as it 
should be and is physically reasonable as explained in [16, 17, 18]. Therefore, this should 
by no means be used as argument against the choice of the energy variable W. 

4.2 The small-Q 2 regime 

To study the dipole formula in the small-Q 2 regime, that is for t — > —00, we rewrite (4.9) 
in the following form 

F 2 (Z, Qle 2t )e- 2t = f°° dt' K (t - t') exp[-(t - t')\ exp(-t')S(£, t') . (4.23) 

J —00 

The function k(t) exp(— t) resembles very roughly a downward step function, see the right 
panel of Figure 11. Furthermore we have 



exp(-t)S(£,t) = \a(r,0 



(4.24) 

r=ro exp(— t) 



and we study here the limit t —> —00, that is, the limit of large r. Let us assume now that 
<r(r, £) is a smooth, only slowly varying, function of r for r — > 00. This is, for instance, the 
case for dipole models where the dipole-proton cross section saturates for r — > 00, 

o~(r,£,) (To for r-> 00. (4.25) 

Now we can apply the same type of argument as we used for the large-Q 2 regime. In 
(4.23) the slowly varying function exp(— t')S(£, t') is, for t — > —00, folded with a function 
resembling a downward step function. Thus, we get an approximate derivative relation 

c_ 00 exp(-t)5(e,t) « -j t [F 2 ^Q 2 e 2t )e- 2t )] , (4.26) 

or, put differently, 



-oc 



4" (r ' " 4^: ( - 2) q2 w [aT{ ^ Q2) + aL ^ Q2)] 



(4.27) 



where 



2 



:= 2tt(k(t) exp(-r))| r ^_ 00 = ^ = 0.093 . (4.28) 

In the right panel of Figure 12 we show for the GBW model the functions of the l.h.s. and 
r.h.s. of (4.26). We see that there is indeed approximate equality of the two functions for 
t < -2. Translated into cross section implies (4.27) to hold for Q 2 < 0.2 GeV 2 . These 
are extremely small Q 2 , and the actual relation is therefore more of academic than of 
phenomenological interest. Nevertheless, we find it useful to discuss its meaning. As we 
will show now, (4.27) gives interesting insight concerning the approximations underlying 
the dipole picture. 
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The result (4.27) is somewhat surprising. It means the following. The dipole model 
with energy variable £ can only work down to very small Q 2 with non-negative cross section 
a(r, £) if or(£, Q 2 ) + 0\l(£, Q 2 ) increases with decreasing Q 2 for small Q 2 . A saturating 
dipole cross section 

a{r, f ) -»• ct / for r -»• oo (4.29) 
implies a logarithmic increase of or + 

^t(£,Q 2 )+^ l (£,Q 2 ) ocln(Qg/Q 2 ) forQ 2 ^0. (4.30) 

Clearly, the behaviour (4.30) is unacceptable physically. For the case £ = W we know that 
gauge invariance requires the behaviour 

<rr(W, Q 2 ) + a L (W, Q 2 ) -> const. for Q 2 -»• . (4.31) 

Indeed, a realistic behaviour for Q 2 — > is 

<tt(W, Q 2 ) + <7 L (W, Q 2 ) = a rp (W) + 0(Q 4 )) , (4.32) 

where a lp {W) is the real-photon-proton cross section and m 2 ff > is a constant of dimen- 
sion mass squared. Inserting (4.32) in (4.27) leads to 

for r — > oo. 

Thus, we have two options. We can either assume that the naive dipole picture with 
the standard perturbative photon wave functions holds down to very small Q 2 values. 
Then, the above arguments force us to give up the saturation hypothesis (4.29) for the 
dipole-proton cross section a(r,W) which should then instead vanish as 1/r 2 for r — > oo, 
see (4.33). We think, however, that the more likely resolution of the above puzzle is that 
the naive dipole picture must be modified for small Q 2 . But we know already from the 
discussions in [16, 17, 18] that the naive dipole picture is expected to break down for 
small Q 2 , where the estimate was Q 2 < 2 GeV 2 . Thus, the discussion above gives further 
evidence for this breakdown. 

In [16] a general analysis of the modifications and additions to the photon wave func- 
tions and the dipole cross section to be expected for small Q 2 was given. We believe that 
these corrections are relevant in the context of the discussion above. On a more phe- 
nomenological level, one can also cure the problem by introducing modifications of the 
photon wave functions, as done for instance in [23], and get then a reconciliation of (4.29) 
and (4.31). The implementation of vector meson dominance at small Q 2 has similar effects 
with regard to the problem pointed out here. 

5. Conclusions 

The dipole model is widely used in the analysis of deep inelastic scattering data, but also 
in other processes. Often one attempts to extract subtle effects from the data with the help 



1 ° f ^ 



(4.33) 
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of the dipole picture, for example the presence or absence of saturation effects. Many of 
these studies use Bjorken-x as the energy variable in the dipole cross section a. However, 
general considerations based on quantum field theory lead to the conclusion that the correct 
energy variable in a is W. In particular, the dipole cross section has to be independent 
of the photon virtuality Q 2 . This issue of the energy variable in a has not been taken 
very seriously so far, since the integral over dipole sizes is thought to be dominated by 
dipoles of a typical size r ~ l/Q- Assuming the strong dominance of a typical dipole 
size one can freely change from one energy variable to another, in particular from W to x 
and vice versa. We have shown here that this simple picture is misleading. In particular, 
we have shown that the corresponding change of energy variable has a large effect on the 
resulting structure function for any given model of the dipole cross section. Numerically, 
the resulting effect is sizeable and can easily exceed the spread among different models for 
<j suggested in the literature. For a reliable study of subtle effects in the data the correct 
choice of energy variable hence appears much more important than has been previously 
assumed. 

Apart from the numerical consequences of choosing energy variables other than W in 
the dipole cross section there is also an important conceptual reason for choosing W. With 
the variables used in (2.5) that formula is an actual factorisation: the variables of the l.h.s., 
Q 2 and W, are separated and occur in two different factors in the integrand on the r.h.s. 
If the second factor in the integrand on the r.h.s. would still contain x and hence Q 2 , the 
dipole formula would not constitute a factorisation. This would spoil the familiar (and 
correct) picture of photon-proton scattering as a two-step process: the photon fluctuating 
into a quantum-mechanical superposition of dipoles of all possible sizes and the subsequent 
dipole-proton scattering. The photon virtuality Q 2 determines the probability distribution 
of dipole sizes, but a dipole of given size does not inherit any information on Q 2 . Its 
subsequent scattering on the proton hence cannot depend on Q 2 . If, as is usually done in 
the dipole picture for DIS, the second step is interpreted as the cross section of asymptotic 
dipole states on the proton, the use of an energy variable involving Q 2 in the dipole cross 
section would violate the rules of quantum field theory. 

In the second part of the paper we have shown that the dipole formula can be writ- 
ten such that the structure function F2 is represented by a convolution of the (suitably 
rescaled) dipole cross section with a known function originating from the photon wave 
functions. In principle, this convolution could even be invertible, albeit with some caveats. 
While a precise determination of the dipole cross section from F2 based on inversion of 
the convolution appears very difficult we have found that the asymptotic behaviour of the 
dipole cross section <r(r, £) for small and large r can be extracted from F2 (even for any 
given choice of energy variable £). In particular, we have found derivative relations between 
F2 and the dipole cross section for both large and very small Q 2 . We expect that these 
findings can be used to construct improved models for the dipole cross section. 

Further, we have been able to explain why in general dipole models using the energy 
variable x are better suited for fitting F2 in a wide range of kinematic parameters than 
models using the correct variable W. Based on the representation of the dipole formula 
as a convolution this can be directly traced back to the sign of the scaling violations that 
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occur for increasing Q 2 at fixed small x. In our opinion that does not disfavour the choice 
of W as energy variable in the dipole cross section. On the contrary, the dipole model 
clearly contains approximations implying already a priori a limited kinematic range for its 
applicability. As shown in [16, 17, 18] the choice of W as energy variable leads to such 
limits which are physically reasonable. 
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A. Integrated photon densities 

In this appendix we derive the asymptotic behaviour of the (over longitudinal momentum 
fraction a) integrated photon densities w^ L (r, Q 2 ), defined in (2.6) and (2.7), at small and 
large distances r. In the main text these limits are needed only for massless quarks. For 
completeness we give them for massive quarks as well. 

At small distances, rQ <C 1 and rm q <C 1, we find from (2.8) and (2.9) 



w^(r,Q 2 ) 
w ( ?(r,Q 2 ) 



>0 N ca C mQl 



2tt 2 



1 



3 (Qr) 2 ' 



2N c a cm Q 2 _ 1 



7T Z 



Q'-\n 2 (Qr). 



(A.l) 
(A.2) 



Note that these leading terms are independent of m q . 

For the limits at large distances we consider massless quarks and massive quarks sep- 
arately. For massless quarks, m q = 0, we find without any approximation 
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(A.4) 



with Meijer's G- function. From this we get at large distances, rQ S> 1: 
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For massive quarks, m q > 0, and large distances, rQ S> 1 and rm q 3> 1, our calculation 
for Wy' and conjecture for to^ give 



w, 



(?) 



(9) 



(r,Q 2 



(r,Q 2 ) 



A^aemQ 2 , n 2 7Tm g exp(— 2m 9 r) 
2vr2 ^ " _ 

2 ^ 



2m„ 



T](m q /Q) 



exp(—2m q r) 



(A.7) 
(A.8) 



with some r-independent function r](m q /Q) for which we find in the case m q < Q the 
numerical value rj(m q /Q) ~ 0.25 independent of m q /Q. The formula (A.8) is based on the 
assumption that for large r the factor multiplying exp(— 2m q r) can be expanded in a series 
in 1/r, and the leading exponent has been determined numerically. 

Figure 13 shows the photon densities as functions of the dipole size r together with 
the asymptotic expressions given above. 
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Figure 13: Integrated photon densities Wjf 1 (r, Q 2 ) (upper plot) and w£ \r, Q 2 ) (lower plot), both 
normalised to a cm Q 2 , as a function of the dipole size r. In addition, the leading terms for small 
and large r are shown as dotted curves. The photon virtuality is fixed to Q 2 = 10 GeV 2 . 
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B. A simplified version of the Donnachie-Dosch model 

In section 3 we have used a simplified version of the Donnachie-Dosch model of [29]. In 
this appendix we briefly describe how the simplified model used in the present paper is 
obtained from the original model. 

We start from eq. (3) of [29] which is the cross section for the scattering of two dipoles 
of sizes R\ and R 2 , 

a Aiv {R u R 2 ) = 0.67 l-((g 2 FF)a 4 f R ± (l - e~^) R 2 (l - e"^) , (B.l) 

with the parameter a = 0.346 fm and the gluon condensate (g 2 FF) taken from lattice 
results, {g 2 FF)a A = 23.77. The cross section for photons or hadrons as external particles 
is obtained by folding the dipole-dipole cross section with the respective wave functions. 
This is given in [29] for the process ab — > cb. For our case of the elastic amplitude 7*p — > j*p 
these wave functions are those of the photon and of the proton. In [29] the latter is 

MR) = 7^r- p exp ("5) 

with Rp = 0.75 fm. Combining formulae (2) and (5) of [29] and relabelling R\ as our r this 
leads to our formulae (2.5)-(2.7) with 

a(r) = J d 2 R 2 \MR2)\ 2 a dip (r,R 2 ). (B.3) 

In addition, an energy dependence of the dipole cross section is introduced by hand 
in [29]. It represents the exchanges of a soft and of a hard pomeron. The soft pomeron 
contributes only if both dipoles are larger than a certain R c , chosen to be R c = 0.22 fm. 
The hard pomeron contributes only if at least one of the two dipoles is smaller than R c . 
This leads to four different integration regions in the integrals over the two dipole sizes, 
according to whether r and/or R 2 are smaller or larger than R c . In our simplified model 
we assume instead that the soft pomeron is the only contribution if r (the size of the qq 
pair originating from the photon) is larger than R c , and the hard pomeron is the only 
contribution if r < R c , independently of the size of R 2 relative to R c . 

Taking into account only these two contributions with their assumed energy behaviour 
we have, similar to eq. (14) of [29], 

fW\ 2es fW\ 2eh 
a(r, W) = 6(r - R c ) a s (r) ^—j + 9(R C - r) a h {r) J (B.4) 

with the parameters given in (3.13). The theta- functions indicate where the two contribu- 
tions are relevant. Here a s is given by the independent factors of eq. (16) of [29], that is 
by (B.3) above. 2 The contribution &h, on the other hand, is obtained from eq. (17) of [29] 

2 Note that there is a factor R1R2 missing in the integrals of equations (16) and (17) in the eprint-version 
of [29]. The journal version contains these factors. We thank H. G. Dosch for clarifying discussions of this 
point. 



(B.2) 
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by dropping the integral over Ri and the photon wave functions. As explained above we 
take into account only the contribution where r < R c , but integrate over all dipole sizes 
in the proton. In eq. (17) of [29] this corresponds to only the second of the three integrals 
there, with the lower limit of the ^-integration set to zero. Accordingly, 



The integral over R2 occurring in (B.3) and (B.5) can be performed numerically. Concen- 
trating on the i?2-dependent factors only we have 



Collecting all factors, we finally arrive at the simplified model (3.12) with the parameters 



Finally, the original model of [29] introduces a rescaling of the dipole cross section by 
the running coupling a s (Q 2 ), see eq. (13) there. That rescaling is relevant in particular at 
large Q 2 . We leave out this factor in order to avoid a Q 2 -dependence of the dipole cross 
section. 
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